B° — B° mixing beyond factorization 
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We present a calculation of the B° — B° mixing matrix element in the framework of QCD sum rules 
for three-point functions. We compute a s corrections to a three-point function at the three-loop level 
in QCD perturbation theory, which allows one to extract the matrix element with next-to-leading 
order (NLO) accuracy. This calculation is imperative for a consistent evaluation of experimentally- 
measured mixing parameters since the coefficient functions of the effective Hamiltonian for B° — B° 
mixing are known at NLO. We find that radiative corrections violate factorization at NLO; this 
violation is under full control and amounts to 10%. 
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The phenomenon of particle- antiparticle mixing, possi- 
ble in systems of neutral mesons of different flavors, is the 
primary source of studies of CP violation (for review, see 
e.g. Q). According to the Cabibbo-Kobayashi-Maskawa 
(CKM) picture, quarks of all three generations must be 
present in a transition for CP violation to occur. His- 
torically, studies of K° — K° mixing provided first essen- 
tial insights into the physics of heavy particles as well 
as tests of general concepts of quantum field theory. For 
a long time it was the only place where the effects of 
CP violation were clearly established (see e.g. 0). Since 
weak couplings of s and d quarks to third generation 
quarks are small, experimental studies of CP violation 
in heavy mesons are considered more promising. While 
recent experimental results for heavy charmed mesons 
D(uc) are encouraging, a full consistent theoretical de- 
scription of this system is still lacking Q ■ These consid- 
erations make the systems of Bd{db) and B s (sb) mesons 
the most promising laboratory for a precision analysis of 
CP violation and mixing both experimentally and theo- 
retically Q. Hereafter we shall consider B d mesons. The 
generalization to B s mesons is straightforward. 

Phenomenologically the system of the B° — B a mesons 
is described by the effective mass operator (M — iT/2)^, 
= {1,2} which in the presence of AB = 2 inter- 
actions acquires non-diagonal terms. The difference be- 
tween the values of the mass eigenstates of B mesons 
Am = Mheavy - M Ught ps 2|Mi 2 | is an important 
observable which is precisely measured to be Am — 
0.489 ± 0.005(stat) ± 0.007(syst) ps' 1 0. With an ad- 
equate theoretical description, it can be used to extract 
top quark CKM parameters. 

In the standard model, the effective low-energy Hamil- 
tonian describing AB — 2 transitions has been computed 
at next-to- leading order (NLO) in QCD perturbation the- 
ory (PT) @ 
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where rjg = 0.55±0.1 J5 — 1.627 in the naive dimen- 
sional regularization (NDR) scheme, So(x t ) is the Inami- 
Lim function |8j, and 0(p) = (&i,7 CT (ii,)(5i,7 -(ii)(/u) is a 
local four-quark operator at the normalization point \x. 
Note that the part of Eq. Q in the second line is re- 
normalization-group (RG) invariant. Mass splitting of 
heavy and light mass eigenstates can then be found to 
be 



Am = 2\(B°\H^= 2 \B°)\ 
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where C = Gj^M^ (V tb *V td ) 2 r,Bm B S {xt)/ (4^ 2 ) . The 
largest uncertainty of about 30% in the theoretical cal- 
culation is introduced by the poorly known hadronic ma- 
trix element A = (B°\O{n)\B ) Q. The evaluation of 
this matrix element is a genuine non-perturbative task, 
which can be approached with several different tech- 
niques. The simplest approach ("factorization") re- 
duces the matrix element A to the product of matrix 
elements measured in leptonic B decays A? = (8/3) 
{B°\bL-fadL\0)(0\bL-f a d L \B o ) = (2/3)_/fm| where the 
decay constant Jb is defined by (0|6l7^c?l|-B (p)) = 
A deviation from the factorization ansatz is 
usually described by the parameter Bb defined as A — 
BbA^; in factorization Bb — 1. There are many ap- 
proaches to evaluate this parameter (and the analogous 
paramete r Bk of K° — K° mixing) available in the liter- 
ature |Ii[Il!]3[i!QBB[!3 

The calculation of the hadronic mixing matrix ele- 
ments using Operator Product Expansion (OPE) and 
QCD sum rule techniques for three-point functions llll 
Iha . Il3j | is very close in spirit to lattice computations jlfij . 
which is a model-independent, first-principles method. 
In the QCD sum rule approach one relies on asymptotic 




FIG. 1: Perturbation theory diagram at LO 
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FIG. 2: Factorizable diagrams at NLO 



expansions of a Green's function while on the lattice the 
function itself can be computed in principle. The sum 
rule techniques also provide a consistent way of taking 
into account perturbative corrections to matrix elements 
which is needed to restore the RG invariance of physical 
observables usually violated in the factorization approx- 
imation (ljj . The calculation of perturbative corrections 
to B° — B° mixing using OPE and sum rule techniques is 
the main subject of this paper. A concrete realization of 
the sum rule method applied here consists of the calcula- 
tion of the moments of the three-point correlation func- 
tion of the interpolating operators of the _B-meson and 
the local operator O(p) responsible for B° — B° transi- 
tions. 

Let us consider the three-point correlation function 

n(pi,p 2 ) = J dxdy(0\TJ s (x)O(0)J B (y)\0y P2X - ipiy . 

The operator Jg = (m& + m^dvy^b is chosen as interpo- 
lating current for the i? -meson and is the b quark 
mass. Note that Jb is RG invariant, Jb = d^(d r y^ L r y^b) 
and (0| Jb(0)\B°(p)) — Jb^b wnere m B is the B- meson 
mass. A dispersive representation of the correlator reads 



n( Pl , P2 ) = n( P lplq 2 



p{s 1 ,s 2 ,q 2 )ds 1 ds 2 
(si -p\){s 2 -p\) 
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where q — p 2 — p\. For the analysis of B° — B° mixing 
this correlator needs to be computed at q = 0, while 
within the sum rule framework q 2 = 0. This particular 
kinematical point is infrared safe for massive quarks. The 
matrix element (B \O{p)\B°) appears in the three-point 
correlator as a contribution of the £?-mesons in the form 
of a double pole 
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where the ellipsis stand for higher resonances and contin- 
uum contributions. The matrix element can be extracted 
by comparing the representations given in Eq. Q and 
the (smeared) theoretical expression of Eq. © obtained 
with an asymptotic expansion based on OPE. Note that 
the analytical calculation of the spectral density itself at 
NLO of PT expansion is beyond present computational 
techniques. Therefore, a practical way of extracting the 
B° — B° matrix element is to analyze the moments of the 
correlation function at pf = p 2 = at the point q 2 = 0. 
One obtains 



i\j\dp\ l dp 2 



p(si,s 2 ,0)dsids 2 

A l A 2 



A theoretical computation of these moments reduces to 
an evaluation of single scale vacuum diagrams (we ne- 
glect the light quark masses). This calculation can be 
done analytically with available tools for the automatic 
computation of multi-loop diagrams. 

The leading contribution to the asymptotic expansion 
is given by the diagram shown in Fig. ^ At leading or- 
der (LO) in QCD perturbation theory the three-point 
function of Eq. (j3Jl completely factorizes H(pi,p 2 ) — 
(8/3)n M (pi)L™(p 2 ) where n^p) is the two-point corre- 
lator 
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The calculation of moments is straightforward since the 
double spectral density p(s±, s 2 ,q 2 ) can be explicitly 
found. Using a dispersive representation of II(p 2 ) 
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one finds the LO double spectral density p LO (si, s 2 , q 2 ) = 
(8/3) (pi ■ P 2)p(si)p(s 2 ) = (4/3)( Sl + s 2 ~ q 2 )p(si)p( S2 ). 
First non-factorizable contributions to Eq. appear at 
NLO. Nevertheless, the factorizable diagrams form an 
important subset of all contributions, as they are inde- 
pendently gauge and RG invariant. Thus, a classification 
of diagrams in terms of their factorizability is a very pow- 
erful technique in the quantitative analysis. 

The NLO factorizable contributions are given by the 
product of two-point correlation functions from Eq. JSJ, 
as shown in Fig.0 Writing Il(p 2 ) = n LO (p 2 )+n NLO (?> 2 ) 
we obtain H^ LO (p l7 p 2 )^(8/3)(p 1 .p 2 )(n LO (pl )U NLO {pl) 
+riNLo(Pi)nLo(P2))- The spectral density of the corre- 
lator LTnloCp 2 ) i s known analytically. This completely 
solves the problem of the NLO analysis in factorization. 
Note that even a NNLO analysis of factorizable diagrams 
is possible as several moments of two-point correlators are 
known analytically. Others can be obtained numerically 
from the approximate spectral density [lflj . 

The NLO analysis of non-factorizable contributions 
within perturbation theory is the main result of this pa- 
per. This analysis amounts to the calculation of a set 
of three-loop diagrams (a typical diagram is presented 
in Fig. I3J ■ These diagrams can be computed using the 
package MATAD for automatic calculation of Feynman 
diagrams . Before applying this package, the combi- 
natorics of disentangling the tensorial structures has to 
be solved and all the diagrams have to be reduced to a 
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FIG. 3: An example of a non-factorizable diagram at NLO 



set of scalar integrals which can be done using the results 
of ref. 2lj. The steps described above were automized 
with the computer algebra system FORM 22] . We shall 
present the details of this calculation elsewhere. 

The local four-quark operator O entering the effective 
Hamiltonian has to be renormalized. We employ dimen- 
sional regularization with an anticommuting 75. The 
renormalization of the operator O reads 



QH = QL 
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with O c = (b L T ^ a t a d L ){b L T^ a t a d L ). The t a are the 
SU C (3) generators and T^ va = (j /J ,^ v j a -j a ^ fl )/2. The 
renormalization of the factorizable contributions reduces 
to that of the 6-quark mass m. We use the quark pole 
mass as a mass parameter of the calculation. 

The expression for the "theoretical" moments reads 



M th (i,j) 



m 2(i+j) 
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where the quantities a^-, b{j and b™? represent LO, NLO 
factorizable and NLO nonfactorizable contributions as 
shown in Figs. Ill.il The NLO nonfactorizable contribu- 
tions 6™^ with i + j < 7 are analytically calculated in 
this paper for the first time. The calculation required 
about 24 hours of computing time on a dual-CPU 2 GHz 
Intel Xeon machine. The calculation of higher moments 
is feasible but requires considerable optimization of the 
code. This work is in progress and will be presented else- 
where. As an example, we give the analytical results for 
the lowest finite moment M t ^(2,2): 
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Here S 2 = ^Cl 2 



= 0.2604..., Cs = C(3), and 



u, = m . For higher moments we present only numeri- 
cal values for 6™/ : 62(2345) = {0.68,1.22,1.44,1.56} and 

b n J 



u 3(34) ~~ {1-96,2.25}. The nonperturbative contribution 
due to the gluon condensate is small. For the standard 
numerical value {g 2 s G 2 ) — 0.5 GeV 4 ,23j], the nonfactor- 
izable contribution due to the gluon condensate amounts 
to about 3% of the perturbative contribution for the high 
moment M(3, 3). For lower moments it is even smaller 



and, therefore, can be safely neglected in the whole anal- 
ysis. 

We use the above theoretical results to analyze sum 
rules and extract the non-perturbative parameter Bb- 

The "phenomenological" side of the sum rules is given 
by the moments which can be inferred from Eq. Q 



M ph (i,j) = -B 
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where the contribution of the i?-meson is displayed ex- 
plicitly. The remaining parts are the contributions due 
to higher resonances and the continuum which are sup- 
pressed due to the mass gap A in the spectrum model. 

For comparison we consider the factorizable approxi- 
mation for both "theoretical" 



<(U) = 
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and "phenomenological" moments, which, by construc- 
tion, are built from the moments of the two-point func- 
tion of Eq. © 
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According to standard QCD sum rule technique, the 
"theoretical" calculation is dual to the "phenomenologi- 
cal" one. Thus, Eq. (fTT7f> should be equivalent (in the sum 
rule sense) to Eq. (JHJl. Also, in factorization, Eq. (|12fl is 
equivalent to Eq. (JTTJ)- Now Eq. JTTJ and Eq. © differ 
only due to non-factorizable corrections. Therefore, the 
difference between Eq. I|12f) and Eq. (|1U|) is because the 
residues differ from their factorizcd values. 

To find the nonfactorizable addition to Bb from the 
sum rules we form ratios of the total and factorizable 
contributions. On the "theoretical" side one finds 



M th (i,j) 



unf 



47T 1 4- 2« hf. 



(13) 



This ratio is mass-independent. On the "phenomenolog- 
ical" side we have 



M ph (i,j) = B B + Rb{z 3 + z l ) + C B z i+i 
M f h l + B/{zJ +z i ) + Cfz i +i 



(14) 



where z — m 2 B /(rn 2 B + A) is a parameter that describes 
the suppression of higher state contributions. A is a 
gap between the squared masses of the B-meson and 
higher states. Rb, Cb, R f and are parameters of 
the model for higher state contributions within the sum 
rule approach. In order to extract the non-factorizable 
contribution to Bb we write Bb = 1 + AB. Similarly, 
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one can parameterize contributions to "phenomenolog- 
ical" moments due to higher B-meson states by writ- 
ing R B = Rf + AR and C B = C* + AC. Clearly, 
AS = AR = AC = in factorization. We obtain 

M ph {i,j) _ 1 | AB + AR(z> + z*) + ACV+J 
Af^,(i,i) l + Rfizi + ^ + Cfzi+i 

Comparing Eqs. I|13(l and (|15(l one sees how the per- 
turbative non-factorizable correction 6™ is "distributed" 
among the phenomenological parameters of the spec- 
trum. We extract AB by a combined fit of several "the- 
oretical" and "phenomenological" moments. The final 
formula for the determination of AB reads 

= AB + AR(z j ~ 2 + z 1 - 2 ) + ACz i+j - 4 (16) 

where AR and AC are free parameters of the fit. We 
take A = 0Am B which corresponds to the duality in- 
terval of 1 GeV in energy scale within the HQET anal- 
ysis of the B meson two-point correlator. We then 
perform least-squares fit to determine AB. Using all 
available theoretical moments we find (AB, AR, AC) = 
a s (m)/(47r) (7.1, -5.0,3.6). We checked the stability of 
the sum rules which lead to a prediction of AB. It can 
be illustrated in the following way. The contribution of 
higher B states is suppressed more strongly for higher 
moments and therefore decreases with increasing order 
of a moment, while the perturbative correction grows. 
The sum of both is (approximately) the same for all mo- 
ments, which leads to a (almost) constant value for AB, 
independent of the particular moment. The calculation 
can be further improved with the evaluation of higher mo- 
ments. The result is sensitive to the parameter z or to the 
magnitude of the mass gap A used in the parametriza- 
tion of the spectrum. Estimating all uncertainties we fi- 
nally find the NLO non-factorizable QCD corrections to 
AB due to perturbative contributions to the sum rules 
to be AB = (6 ± l)a s (m)/(47r). For m = 4.8 GeV, 
a s (m) = 0.2 H |3 it leads to AB = 0.095 w 0.1. 
It is known that nonperturbative corrections (such as 
the ones due to the quark-gluon condensate) to the pa- 
rameter B B are negative, AB nonPT {m) = -0.05 [l^ . 
Combining this result with the present analysis we find 
B B (m) — 1 + O.Ipt — 0.05„ o „pt showing the excellent 
numerical validity of the factorization approximation at 
the scale \i — m. 

In conclusion, we have evaluated the B° — B° mixing 
matrix element in the framework of QCD sum rules for 
three-point functions at NLO in perturbative QCD. The 
effect of radiative corrections on B B is under complete 
control and amounts to approximately +10%. We have 
also shown that perturbative QCD correction to AB for 
the moments considered in our analysis completely dom- 
inates the correction due to the gluon condensate. 
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